The coil-globule transition of polymers of long rigid monomers connected by flexible 

spacers 



Richard P. Sear 

FOM Institute for Atomic and Molecular Physics, Kruislaan 407, 1098 SJ Amsterdam, The Netherlands 

email: sear@amolf.nl 

A simple model of a polymer with long rigid segments which interact via excluded volume 
repulsions and short ranged attractions is proposed. The coil-globule transition of this model 
CT^N ' polymer is strongly first order, the globule is crystalline and the coil which coexists with the globule 

is swollen. A virial expansion truncated at low order is shown to provide a very poor approximation 
to the free energy and so a cell theory is used to calculate the free energy of the globule. 
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I. INTRODUCTION 



Typically, in a dilute polymer solution, each polymer chain is spread over a very large volume.EJ The polymer chain 
forms a fractal and if a box is drawn so as to enclose the chain, almost all of the volume within the box is occupied 
by solvent, not by the monomers of the polymer. This is the coil state of a polymer chain. The coil may be ideal, in 
which case its radius of gyration scales with the number of its monomers N as N 1 / 2 , or it jnay be swollen, in which 
casa-its radius of gyration scales as N 3 ^ 5 . The value 3/5 of the exponent is the Flory value,El the true value is slightly 
less.til However, there is a second state of the polymer: the globular state, in which the radius of gyration scales as 
iV 1 / 3 . As the radius of gyration scales with N 1 ^ 3 , the volume scales linearly with N, the density is then finite (except 
at a second order transition to the coil state). The volume occupied by a coil increases more rapidly than linearly 
with N and so as N becomes very large the density tends to zero. The state of a polymer is in turn determined by 
the interactions between different parts of the polymer. If they repel each other then the polymer is the swollen coil 
state, if they attract each other strongly enough then the polymer is in the globular state, whilst if the attractive and 
repulsive parts of the interactions balance the polymer is ideal apart from logarithmic corrections.lirEj By decreasing 
the temperature and so increasing the effect of the attractive forces the polymer may be transformed from a coil to a 
globule. A considerable amount of theoretical work has been devoted to attempting to understand and describe the 
£^ coil-globule transition. However, this has concentratedj-almost exclusively on coil-globule transitions driven by the 
polymer's second virial coefficient becoming negative.mJ In the coil state, the interaction part of the free energy is 
£S) ' dominated by interactions between only two parts of the polymer, which are accounted for by the polymer chain's 
second virial coefficient (which is general different from the second virial coefficient of the monomers). So, as the 
second virial coefficient decreases to zero the interactions effectively vanish and the coil becomes (almost) ideal. Then 
on reducing the temperature further the second virial coefficient becomes negative and the coil collapses into a globule. 
The coil first becomes ideal before collapsing only when the collapse is driven by the change of sign of the polymer's 
second virial coefficient. Here, we show that when the monomers are long and rigid the coil-globule transition is not 
driven by the second virial coefficient becoming negative. So, the polymer is never ideal. When the temperature is 
reduced the polymer transforms-via a strongly first order transition from a swollen coil to a dense crystalline globule. 
■ The use of a virial expansiont£l in order to describe the contribution of monomer-monomer interactions to the free 
energy of the polymer is reasonable for spherical monomers with not too short ranged attractive interactions. However, 
O ■ if the monomers are long and rigid an additional length scale, the length of the monomer L, is present. Unless the 
range of the attraction, is as long or longer than this length L then the attractive forces between the monomers are 
highly anis^ropic.EJil3 Then, their effect on the third and higher virial coefficients is larger than on the second r*4ijial 
coefficient .113 Long rigid monomers behave in a similar way to hard spheres with a short ranged attraction.E-3~li3. 
rN \ StclOEl has shown that, as the range of attraction between the spheres tends to zero, the virial expansion becomes 
pathological at temperatures such that the attractive interactions make a nonnegligible contribution to the second 
virial coefficient. Here, we show that the same occurs for long rods as the range of attraction goes to zero. Unless 
the renormalisation requirecn on going from the virial coefficients of the monomers to the effective virial coefficients 
of a polymer completely changes the behaviour of the virial series, the coil-globule transition cannot be driven by the 
renormalised second virial coefficient becoming negative. i_. _ . 

The DNA double helix is a rather rigid polymer, its persistence lengthEj is ~ 50nmt3 r^hich is 25 times its diameter 
of ~ 2nm. The interaction between DNA helices has a much shorter range than SOnmElJE 3 ] and so we would expect 
the behaviour of a DNA molecule of length 3> 50nm, to be qualitatively similar to that of our mod^L-polymer chain. 
This expectation is borne out, the coil-globule transition of DNA, as for our model, is first order .1 
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II. THE MODEL POLYMER CHAIN 



The polymer is a linear chain of N cylindrical monomers joined with (N — 1) flexible spacers, see Fig. (yj). A spacer 
is flexible enough to allow the two cylinders which it joins to adopt any relative orientation with equal probability 
and to lie side-by-side, however, its maximum length is much less that the length L of a cylindrical monomer. We 
will always consider N to be infinite and so the coil-globule transition is a phase transition. Real polymers are, of 
course, always finite but we do not consider the effect of finite size here. 

The cylindrical monomers are of length L, and have a hard core of diameter D; when two cylinders overlap the 
energy of interaction is positive and infinite. The attractive part of the potential is sticky in the sense that its range 
is close to zero, it is only nonzero between a pair of cylinders if they almost touch for all of their length, for this to be 
true they must be almost parallel and their centres of mass must be almost side-by-side, i.e., only a little more than 
D apart. We will define what we mean by almost below when we define the potential. When the attractive part of 
the potential is nonzero it equals — e. The interaction u between two cylinders, f and 2, is a function of the separation 
ri2 of the centres of mass of the cylinders, and of the angle 712 between the centrelines of the two cylinders. Then, 
w(ri2,7i2) is defined by 

if D < < (1 + 5)D, \zi 2 \<SD, j 12 < 8 u = -e 

else if the hard cores overlap u = 00 (f ) 

otherwise u = 0, 

where the vector r 12 = (r\ 2 , Z12, 4>i 2 ) in cylindrical polar coordinates. The origin of the coordinates is fixed on the 
centre of mass of cylinder f and their z-axis is along its axis. The parameter 8-^1 defines how far apart the cylinders 
can be and still interact via the attractive part of their interactions. T he , . li mit 8 — > defines a zero range potential, 
a generalisation to a cylinder of the sticky sphere potential of Baxter.E-§E3 The dependences of the attractive part 
of the potential u on 712 and on the components of ri2 are independent; this is rather artificial but allows a simple 
evaluation of the second virial coefficient and as 8 is near zero the details of the attractive interaction make little 
difference to the behaviour of the polymer. 

In Eq. (|I|), 8 and e are just adjustable parameters. However, for an attractive interaction-of range ~ D between 
rods of length L the energy of interaction of two parallel rods scales linearly with L/D.li3Eil This produces a very 
deep well, or a large e in terms of our model. The energy continues to scale linearly with L/D as long as the angle 
between the two rods is < D/L. Thus, the L/D — > 00 limit of the potential models considered by van der Schoot and 
OdijkEa is very similar to the e — * 00, 8 — ► and eS -^constant limit of Eq. ([!]). 

III. THE COIL STATE 

The density of a coil tends to zero as N tends to infinity. Thus, the only interactions between distant parts of the 
chain that count are pair interactions; the density i$rSO low that the probability of three different parts of the chain 
coming close enough to interact is vanishingly smahmj By distant parts of the chain we mean monomers which are not 
close together along the chain, i.e., are separated along the chain by many other monomers. As only pair interactions 
are important the excess free energy depends on only one parameter, an effective second virial coefficient B 2 . This 
is not the same as the second virial coefficient between monomers, B 2 . The effective coefficient B 2 is obtained by 
renormalising B 2 \ this 'mixes in' some of the higher virial coefficients in with B 2 to produce B 2 . So, we start by 
determining B 2 and then go on to estimate the higher virial coefficients. 

In the-iimit of the temperature|-tfjnding to infinity, the polymer chain is just a chain of hard cylinders and flexible 
spacers.c3 Then, B 2 = (7r/4)L 2 Z?Jl3'EZI The higher virial coefficients are all very small if L/D I.lScH Thus interac- 
tions between pairs of the monomers are close to being independent, and so if two distant (along the chain) monomers 
are interacting the probability of any of the monomers which are adjacent to either of these monomers interacting 
with either of the original pair of monomers is very low. As it is these interactions which renormalise B 2l we conclude 
that the renormalisation required to derive B 2 from B 2 is weak and so B 2 — B 2 . The pair interactions cause the coil 
to swell so that its radius of gyration scales as A 3 / 5 , not as N 1 / 2 as an ideal coil does. However, they contribute an 
amount much less than T to the free energy per monomer. T is the temperature in energy units. This can be seen if 
we examine the interactions between distant and between adjacent monomers. The interaction between distant pairs 
of monomers, i.e., pairs of monomers which, are not neighbours or next-nearest neighbours etc. along the chain, is 
characterised by the Fixman parameterai3'E3 z ~ B 2 / L 3 N 1 / 2 . If we divide this by N to obtain the contribution per 
monomer we see that the contribution tends to zero as N tends to infinity. As for the interaction between adjacent 
monomers; it is easy to see that a monomer only very weakly restricts the space available to an adjacent monomer. 
Thus, the contribution to the free energy from interactions between monomers is small. Then the free energy of our 
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chain is to a good approximation the free energy of an ideal chain. If the flexible spacers allow each monomer to move 
freely in a phase space of volume D 3 then the free energy per monomer of the coil a is 

a(T) = -T\nD 3 . (2) 

Of course there is also a term from the momenta but we neglect this as it is the same in all phases. Fixing the volume 
of phase space to be D 3 is consistent with short tethers of length not much larger than D. 

In order to consider the effect of attractive interactions we start by studying the second and third virial coefficients 
of sticky cylinders. In fact the contribution of pairwise interactions to the free energy of a coil involves not B 2 but a 
renormalised second virial coefficient B^a The difference betwee n t he two will be considered after the unrenormalised 
coefficients have been calculated. The second virial coefficient iaija 

B 2 {T) = J (exp[-u(r 12 , 7l2 )/T] - 1) dv 12 d^dn 2 , (3) 

where d£li = sin 8 idQidtfii. The diagrammatic representation of B 2 is given by Fig. (||a). See Ref. [lO]for an introduction 
to the theory of diagrams. The integral in Eq. (|^) is straightforward. For 7 > 5 the interaction is just that between 
two hard rods, whose B 2 is well known.EBEZI The contribution of the attractive part of the interaction is evaluated by 
aligning cylinder 1 with the z-axis, then 712 = 9 2 , and as the integration is restricted to small angles sin #2 — 2 . So, 

j-D+SD pSD *8 

B 2 (T) = jL 2 D--exp[e/T} 2vr 12 dr 12 dz 12 2 d6 2 , (4) 

44 Jd J-sd Jo 

where cylindrical polar coordinates have been used for ri 2 . We have assumed that L/D 1 and e/T ^> 1. 
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B 2 (T) = '-L'D - , -D 3 5 i exv[e/T]. (5) 

This expression for B 2 may be compared with those derived by van der Schoot and OdijkEi for more realistic potentials. 
If in Eq. (||), 5 is replaced by D/L and e multiplied by L/D. then it becomes essentially the same as Eq. (3.7) of Ref. 

H 

The twelfth virial coefficient of sticky spheres diverges to minus infinity at alLtemperatures at which the attractive 
interactions make a nonnegligible contribution to the second virial coefficient £3 Here, we show that the third virial 
coefficient -B3 of sticky cylinders diverges to minus infinity at all temperatures for which the attractive interactions 
make a nonnegligible contribution to B 2 , in the limit that S tends to zero. B3 is defined byE3 

B 3 {T) =-— i-3 f Yl (expI-ufry^ijO/ri-lJdriadrisdni^dns, (6) 

•* zj=12, 13,23 

where r 13 is the vector between the centres of mass of rods 1 and 3, and where 713 and 723 are the angles between 
rods 1 and 3, and 2 and 3, respectively. The diagrammatic representation of B3 is given by Fig. (||b). Only the 
scaling of B3 with e and 6 is required, it is 

B 3 (T) - L 3 D 3 -D 6 5 9 exp[3e/T}. (7) 

The first term is from the hard core.@Ei There are three interactions between the three rods, this gives the factor of 
3 in the exponential of the second term of Eq. (0) . The integration over the interaction of the first rod with each of 
the the second and third rods gives a factor of <5 4 , then as the rods are already restricted to be nearly parallel and the 
second and third rods are already restricted to being close to the first rod the integration over the interaction between 
the second and third rods gives only a factor of S. In the limit of 5 — > the second term of B 2 , Eq. (^), is finite if e 
diverges as e — > In S~ 4 . Then the attractive interactions contribute a finite amount to the value of B 2 . Putting this 
dependence of e into Eq. (Q) results in the attractive term in B3 tending to minus infinity as S — > 0. The divergence 
of the twelfth virial coefficient of sticky hard spheres is suppressed by the presence of even very weak polydispersity 
in the diameter of the spheres. In contrast, polydispersity in the diameter D of sticky rods has little effect: B3 still 
diverges. 

Although we cannot consider all virial coefficients we ace able to find a trend in their behaviour. The fourth virial 
coefficient B4 is the sum of the 3 diagrams of Fig. (^c).l!3 The most highly connected of these (with 6 bonds) does 
not correspond to a realisable cluster;oo it is not possible for all of the 4 rods to interact simultaneously with all of 
the other 3 rods via the attractive part of the interaction. Therefore, it is the diagram with 5 bonds (the middle one 
of Fig. (||c)) which dominates at low temperature. Five bonds means a factor of exp[5e/T]. The region of integration 



3 



over which all 5 bonds are simultaneously present scales as <5 14 . This result may be derived if we start from the B3 
integral and add an extra rod, which interacts with 2 of the existing rods. The integral over the 3 rods in B3 yields 
a factor of <5 9 , then restricting the orientation of the fourth rod introduces an additional factor of S 2 . The centre of 
mass of the fourth rod must be restricted in all 3 directions in order that it interact with 2 other rods; this gives a 
factor of S 3 . So, we can find the order of the diagrams in the nth virial coefficient with respect to exp[e/T] and S by 
starting from diagrams in the (n— l)th coefficient and adding another rod such that the cluster is still realisable, i.e., 
that all the rods can simultaneously interact via all the bonds in the diagram. However, this quickly becomes tedious 
as n increases. We have found that the part of B4 due to the attractive interactions scales as exp[5e/T](5 14 , and so 
becomes non-zero at temperatures such that both B2 and B3 are still at their T = 00 limits. The trend seems clear, 
B n starts to differ from its T — 00 limit at a temperature which increases with n. 

So far, we have just studied the virial coefficients of sticky cylinders. As the third virial coefficient diverges at 
all temperatures such that the attractive-interactions influence B2, a fluid of sticky cylinders is unstable at these 
temperatures, just as sticky spheres are.EfJ However, the main object of our study is a polymer chain consisting of 
an sequence of N of these sticky cylinders. In a polymer coil the effective second virial coefficient is not that of the 
monomers B2 but is B2 renormalised by the higher coefficients, B^n The temperature at which the coil is ideal is the 
temperature at which B 2 is zero. At low temperatures the higher virial coefficients are much larger than B2 and so 
it is reasonable to expect the renormalisation to produce a renormalised B\ which differs substantially from the B2 
of Eq. <^). Indeed, there is no reason why the scaling of B\ and B2 with 5 and e should be the same. Unfortunately 
we do not know how to carry out this renormalisation. In conclusion, at high temperatures we can approximate B 2 
by B2 and in the coil state B 2 is all we require. The assumption that ~ B2 near the T = 00 limit is reasonable 
because then for L/D 1 the higher virial coefficients are all very small. At low temperatures we have little idea 
of the size of B 2 and so we are unable to construct a virial expansion of the free energy of a coil. The crossover 
temperature between these two regimes will be estimated in Section V. 



IV. THE GLOBULAR STATE 



A globule, because it is denser than a coil, has a lower entropy. This is compensated by its lower energy due to 
attractive interactions between the monomers of the polymer. Our rodlike monomers only attract each other when 
they are side-by-side and nearly parallel. Given this, it is easy to see that a monomer has the lowest energy when it 
is closely surrounded by 6 other monomers, all of which are nearly parallel and have their centres of mass almost in 
a plane. Then the monomer interacts strongly with 6 other monomers, the maximum number possible. The lowest 
energy state of the polymer is then just a state which allows every monomer to be surrounded in this way by 6 others. 
This occurs in a solid formed from layers of cylinders arranged in 2-dimensional hexagonal lattices. It is not surprising 
that the lowest energy is achieved in a solid but it should be noted that our monomers have to lie side-by-side to 
attract each other and so the energy of any isotropic phase will be much higher than that of a phase with orientational 
ordering. In addition the orientational ordering must be very pronounced, the angle the monomers make with the 
director should not be more than ~ 5. The restriction that the centres of mass must be within (I + 5)D in order for 
them to attract each other forces the density to be very high, roughly a fraction 1/(1 + S) 2 of the density of close 
packing — the maximum possible density. A globule with a volume fraction of, say, 1% would necessarily have a 
much higher energy due to the very short range of the attrac ti o n s. 

As we expect a globule of such high density, a cell theoryli§E!J is used to estimate its free energy. We consider a 
very dense layer of cylinders with their centres of mass on a 2-dimensional triangular lattice; the lattice constant is 
(1 + 5/2)D. As the interlayer interactions are very weak in our model we do not consider the interactions between 
monomers in different layers. Fig. (Jlfb) is a schematic of part of such a layer. In a cell theory calculation, the nearest 
neighbours of a particle are fixed at their places in a perfect lattice and the partition function q\ of a particle trapped 
in the cell formed by them evaluated. So, the centres of mass of the six neighbours of the particle are fixed at positions 
(f + S/2)D away from the centre of the cell of the particle and in the iry-plane of our coordinate system and their 
centrelines are parallel to its z-axis. The cylinder in the cell has 5 degrees of freedom, 3 translational and 2 rotational. 
As our model is defined so that the orientational and translational degrees of freedom are decoupled, estimating the 
volume of phase space available to the rod is straightforward. The orientational degrees of freedom give a factor of 
(l/2)<5 2 , as for B2. We approximate the volume available to the centre of mass of the rod by a hexagonal prism with 
its axis along the z axis. The height of the prism is 28D and the distance from its centre to any of its points is 
(l/2)8D; its volume is then 0.75-y/3(&D) 3 . So, 

qi ~^-D 3 5 5 exp[3e/T] (8) 
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and the free energy per monomer a is 



a(T) = —T In q\ ~ — 3e — T In 



3\/3 



L> 3 <5 5 



(9) 



The integral over the momenta coordinates has been neglected; it is the same in all phases and so has no effect on the 
phase behaviour. In taking a(T) to be the free energy per monomer of our polymer we have neglected the entropy 
associated with the number of ways in which the monomers can be connected by the spacers; this is of order unity. 



V. THE COIL-GLOBULE PHASE TRANSITION 



In Section II, we showed only that the when B 2 becomes negative the higher coefficients are already negative, and 
infinite. Now that we have an estimate of the free energy of the globule we can determine the stability of the coil 
with respect to the globule when B 2 becomes negative. In the limit of S — > oo the second term of B 2l Eq. (||), is 
finite if e diverges as e — > In J -4 . Putting this dependence of e into Eq. (Q) for the free energy of a globule, results in 
the free energy of the globule tending to — oo. Thus, for cylinders with a very short ranged, sticky, attraction the coil 
is unstable with respect to the globule when the (unrenormalised) interaction between a pair of monomers becomes 
attractive. The coil is only stable at temperatures such that B 2 equals its T = oo limit. 

Now, as the free energy per particle of the coil state is constant (at high temperature) and given by Eq. (j^), we 
can estimate the temperature of the coil-globule transition T cg by equating the free energy of the globule state, Eq. 
@, to that in the coil state, Eq. (Q), 

Tcg ~ ~ 3 (10) 



e In ^<5 5 

Above T cg the coil is more stable than the globule, below T cg the reverse is true. At T cg , the coil collapses via a 
strongly first order transition to form a dense globule. This is qualitatively different from the situation when the 
polymer is flexible, there the radius of gyration changes continuously as the temperature is varied.Eruu A posteriori 
justification for the correctness of Eq. ( |l0| ) is provided by the fact that it predicts a T cg at which both B 2 , Eq.(||), 
and B3, Eq. (Q), are equal to their T = 00 limits. So, unless B 2 at T cg contains large contributions from higher order 
virial coefficients the coil is close to its T = 00 limit at T cg and so the free energy per monomer is close to the value 
of Eq. (Q). Physically, large contributions from high order virial coefficients corresponds to clustering, i.e., significant 
numbers of monomers in the coil would be part of clusters of 4 or more monomers.E3 At equilibrium, there seems no 
reason for small clusters to be favoured, if the attractive interactions are strong enough to bind a few monomers into 
a cluster then they are strong enough to form a macroscopic cluster — the globular phase. 

The temperature T cg is an approximation to the, temperature at which a fluid of the (unpolymerised) monomers 
becomes unstable with respect to a dense solid.EilO It is an overestimate as free monomers possess translational 
entropy in the fluid phase which stabilises the fluid. On the basis of our analysis here we expect a fluid of long rodlikc 
objects with short ranged attractions to undergo a strongly first order transition to a very dense phase, when the 
attractions are made sufficiently strong. This appears to happen for short, ~ 50nm, lengths of DNAEJ Wissenburg 
et al.EH have observed that, on increasing the concentration of the DNA, it suddenly forms dense aggregates. 



VI. CONCLUSION 



The polymer never forms an ideal coil, at temperatures above T cg it exists as a swollen coil and below T cg it forms 
a crystalline globule. The coil-globule transition is strongly first order. A cell theory was used to calculate the free 
energy of the globular state and so to find the coil-globule transition temperature T cg ; the virial expansion for the 
free energy of our polymer havmg iieeii shown to be useless at low temperatures. The coil-globule transition of a 
flexible polymer is second order J3lI'l3Ej the radius of gyration varies continuously with temperature. This globule is 
liquid-like, there is no orientational or translational ordering. Clearly, as the stiffness of the monomers is increased 
there is a chaBger-in the character of the coil-globule transition; this change has been observed in simulations of 
lattice,models.a ti3 Experimentally, data is only available for highly flexible polymers, such as polystyreneHoE3 and 
DNA.I§lI, Y}i4ffe^ s scmiflexible The data for flexible chains seems consistent with theoretical results-,based. on virial 
expansionsoO'Q'Ll and with simulations of flexible chains on a lattice.E3c3 However, DNA forms dense&BEjO globules 
in which the DNA helices are hexagonally ordered. We have studied, using theory, a molecular model which also forms 
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dense ordered globules. The major difference between our model and DNA is that DNA is semiflexible, which means 
that the helix bends continuously, whereas our chain consists of completely rigid monomers, which cannot bend at all, 
and spacers, at which the chain can bend freely. We hope to be able to incorporate the effect of continuous flexibility 
in future work. 
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to thank The Royal Society for the award of a fellowship and the FOM institute AMOLF for its hospitality. The 
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FIG. 1. Part of a polymer chain consisting of N rigid cylindrical monomers connected by flexible spacers; a) shows the 
chain in the coil state and b) in the globule state. The cylindrical monomers are shown with a diameter smaller than D in 
order to be able to see the hexagonal packing in the globule state. The curved black lines represent the spacers. 



FIG. 2. The diagrams for the second, a), third, b), and fourth, c), virial coefficients. The black circles represent the rods 
and the lines represent Mayer f functions, which are exp[— u/T] — 1. A diagram with n circles is an integral over the coordinates 
of n rods, interacting by the Mayer f functions shown in the diagram, divided by a symmetry number. 
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